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The time variation of the elastic scattering rate of solar neutrinos with electrons in Super- 
Kamiokande-I was fit to the variations expected from active two-neutrino oscillations. The best 
fit in the Large Mixing Angle solution has a mixing angle of tan^ 6 — 0.55 and a mass squared dif- 
ference of Am^ = 6.3 X 10~^eV^ between the two neutrino mass eigenstates. The fitted day/night 
asymmetry of — 1.8± 1.6(stat)^}'2(syst)% has improved statistical precision over previous measure- 
ments and agrees well with the expected asymmetry of -2.1%. 
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PACS numbers: 14.60.Pq,26.65.+t,96.40.Tv,95.85.Ry 

The combined analysis of all solar neutrino experi- 
ments 0, Q gives firm evidence for neutrino oscillations. 
All data are well described using just two neutrino mass 
eigenstates and imply a mass squared difference between 
Am2 = 3 X IQ-^eY^ and Am^ = 1.9 x IQ-^'eV^ and a 
mixing angle between tan^ 9 = 0.25 and tan^ 6 = 0.65 
This region of parameter space is referred to as the Large 
Mixing Angle solution (LMA) . The rate and spectrum of 
reactor anti-neutrino interactions in the KamLAND ex- 
periment |J] are also well reproduced for these mixing 
angles and some of these Arr? . Over the Ato^ range 
of the LMA, solar *B neutrinos are ^100% resonantly 
converted into the second mass eigenstate by the large 
matter density inside the sun (5|. Therefore, the sur- 
vival probability into Vf, is « sin^ Q. However, due to 
the presence of the earth's matter density, the oscillation 
probability at an experimental site on earth into dif- 
fers from sin'^ Q during the night. Since our experiment 
is primarily sensitive to I'e's, this induces an apparent 
dependence of the measured neutrino interaction rate on 
the solar zenith angle (often a regeneration of Ve 's during 
the night). We employ a maximum likelihood fit to the 
expected solar zenith angle dependence on the neutrino 
interaction rate. Herein, the statistical uncertainty is re- 
duced by 25% compared to our previous measurement 
of the day /night asymmetry 0,0 which consists of two 
flux measurements in two separate data samples (day and 
night). It would require almost three more years of run- 
ning time to obtain a similar uncertainty reduction with 
the previous method. 

Super-Kamiokande (SK) is a 50,000 ton water 
Cherenkov detector described in detail elsewhere . SK 
measures the energy, direction, and time of the recoil 
electron from elastic scattering of solar neutrinos with 
electrons by detection of the emitted Cherenkov light. 
Super-Kamiokande started taking data in April, 1996. 
In this report, we analyze the full SK-I low energy data 
set consisting of 1496 live days (May 31^*, 1996 through 
July 15*l^ 2000). 

The solar neutrino interactions are separated from 
background events by taking advantage of the strong for- 
ward peak of the elastic scattering cross section. The ar- 
rival time of each solar neutrino candidate defines a solar 
direction. Using this direction, we calculate the angle 
^sun between the reconstructed recoil electron direction 
and the solar direction. The data sample is divided into 
A'bi,, = 21 energy bins: 18 energy bins of 0.5 MeV between 
5 and 14 MeV, two energy bins of 1 MeV between 14 and 
16 MeV, and one bin between 16 and 20 MeV. We use 
two types of probability density functions: p(cos ^aun, K) 
describes the angular shape expected for solar i/g's of en- 
ergy E (signal events) and Ui(cos6'sun) is the background 
shape in energy bin i. Each of the rii events in energy 
bin i is assigned the background factor — Ui{cos9iK) 
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FIG. 1: Allowed Regions for the Eccentricity-Induced Sea- 
sonal Solar Neutrino Flux Variation at 68% (dark gray), 95% 
(gray), and 99.73% (light gray) C.L.. 

and the signal factor Si^ — p{cos6iK., E^,). The likehhood 

is maximized with respect to the signal S and the 21 
backgrounds Bi. MC^ is the number of events expected 
in energy bin i using the flux and spectrum of *B and 
hep neutrinos. 

A simple determination of the day/night asymmetry 
is obtained by dividing the data sample into day and 
night and fit cos6^^„ to each sample separately. From 
the obtained day (D) and night (N) rates we calcu- 
late the asymmetry Aj^^ = {D - iV)/(0.5(£> + N)) = 
— 2.1 ±2.0(stat)^]^'2(syst)% which is consistent with zero. 
To take into account time variations in the likelihood fit, 
the signal factors are modified to Si^ = p{cos9iK, Ek.) x 
Zi{a,tf^) where is the event time and a is an ampli- 
tude scaling factor. As a simple example, we measure the 
earth's orbital eccentricity. Since the neutrino fiux is pro- 
portional to the inverse square of the distance between 
sun and earth, the eccentricity induces a seasonal time 
variation. Below 6.5 MeV the background rates can fiuc- 
tuate at time scales of several weeks or longer mainly due 
to changes in the radon contamination in water. There- 
fore, these energy bins are excluded from the eccentricity 
analysis by setting Zi to 1. To measure both the phase 
and amplitude of the variation, both the eccentricity and 
the perihelion is varied around the known values (1.7% 
and ^ January 3rd). Figure shows the allowed ranges 
of parameters at 68%, 95%, and 99.73% C.L.. We mea- 
sure the perihelion shift to be 13 ± 17 days (consistent 
with zero) and the amplitude of the neutrino flux vari- 
ation to be 1.51 ± 0.43 (consistent with one) times the 
amplitude expected from 1.7% eccentricity. 

To study neutrino oscillation-induced time variations, 
we correct for this seasonal effect (with the nominal 
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FIG. 2; LMA Solar Zenith Angle Variation Shapes. The pre- 
dictions are for Am^ = 6.3 x 10~^eV^ in the energy bins 16 
to 20 MeV (top), 12.5 to 13 MeV, 10 to 10.5 MeV, 7.5 to 8 
MeV, and 5-5.5 MeV (bottom). 



perihelion and eccentricity). After that, the additional 
seasonal amplitude variation is 0.48 ± 0.43 times the 
eccentricity-induced variation which is consistent with 
zero. We search for solar zenith angle variations (em- 
ploying the solar zenith angle as the time variable) and 
additional seasonal variation due to the oscillation phase 
(using the distance between sun and earth). In each bin 
i we calculate the rate ri{t) (oscillated Monte Carlo). 
From this rate and the live-time distribution the av- 
erage (r°^), day, and night rates and subsequently the 
day/night asymmetry Ai are computed. Using the day 
(night) live-times Lo (ijv) and the live-time asymme- 
try Ldn = {Ld — LN)l{Q.b{LD + Lf^j), the effective 
asymmetry parameter = 0.25AiL£)jv is computed and 
Zi(a,t) is defined as 



l + a{{l + a,)n{t)/rr -I) 
1 + a X ai 



so that r-(a,t) = Zi{a,t) x r"" has the same average to- 
tal rate rf, but the day/night asymmetry is Ai x a. 
In particular, r^(0,t) = r"" is independent of t and 
r^(l,t) = ri{t). Figure [3 shows the expected solar zenith 
angle variation shapes Zi(l, cos 9z) in five different energy 
bins using an LMA solution and the density model of the 
earth 0. 

The resulting likelihood function is maximized with 
respect to signal S, the backgrounds Bi, and the asym- 
metry scaling parameter a. For the best-fit LMA os- 
cillation parameters (which will be described later) we 
find a — 0.86 ± 0.77 which corresponds to the day/night 
asymmetry 

^DN = -f.8± 1.6(stat);i;2(syst)% 

where —2.1% is expected for these parameters. The sta- 
tistical uncertainty is reduced by 25% with this likelihood 
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FIG. 3: LMA Spectrum (top) and D/N Asymmetry (bottom). 
The predictions (solid lines) are for tan^ 9 = 0.55 and Am^ = 
6.3 X IQ-'^eY^ with <^8s = 0.96xStandard Solar Model 
and (phep ~ 3.6 X Standard Solar Model. Each energy bin is fit 
independently to the rate (top) and the day/night asymmetry 
(bottom). The gray bands are the ±1(t ranges corresponding 
to the fitted value over the entire range 5-20 MeV: A = — 1.8± 
1.6%. 
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FIG. 4: SK Day/Night Asymmetry as a Function of Am^. 
The solid line is expected from two-neutrino oscillations, the 
band (±lcr) results from the fit to the SK data. The mixing 
angle tan^ 9 = 0.55 is used. Overlaid are the allowed ranges 
in Am^ (cross-hatched bands) from the KamLAND experi- 
ment Q. The second band (LMA-I) is favored. 



analysis; however, the resulting day/night asymmetry is 
still consistent with zero. Figure 13 shows the fitted rate 
(top), as well as the day/night asymmetry (bottom) for 
each energy bin separately. The oscillation expectations 
are indicated by the solid lines. The asymmetry fit value 
and uncertainty depends on the solar zenith angle varia- 
tion shapes Zi(l,t) which in turn depend on the oscilla- 
tion parameters. Figure 01 shows the expected day/night 
asymmetry and fit results for each Am^ in the LMA re- 
gion with the best-fit mixing angle tan^ = 0.55. The 
expected day/night asymmetry and the ±lcr band of the 
fit overlap between 5 — 12 x 10~^eV^. 
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FIG. 5: Excluded (SK spectrum and time variation; dark 
gray) and Allowed (SK spectrum, rate, and time variation; 
light gray) at 95% C.L.. Overlaid are the areas excluded just 
by the day /night and seasonal variation (hatched regions in- 
side thick black lines). The graphs at the top (and right) show 



FIG. 6: Allowed area at 95% C.L from the combination of SK 
and SNO (gray) and all solar data (dark gray). The graphs 
at the top (and right) show the difference as a function of 
tan^e (Am^) only: the dashed line is the SK/SNO fit, the 
solid line includes all solar data. The best fit to all solar data 



the X difference as a function of tan 6 (Am ) alone where is tan 6 — 0.42 and Am — 6.0 x 10 ''eV 
the Am^ {tan'^ 9) is chosen to minimize x^- 



To constrain neutrino oscillation using the SK rate 
time variations, the likelihood difference Alog£ = 
log£(a = 1) — log£(a = 0) between the expected time 
variation and no time variation is computed. Below 
Am^ = l.SxlO^^eV^, the day/night variation is replaced 
by an additional seasonal variation due to the oscillation 
phase. As for the eccentricity-induced variation, the en- 
ergy bins below 6.5 MeV are excluded from the seasonal 
variation because of the slow time variation of the back- 
ground. However, since the effect of that variation on 
the day/night asymmetry was carefully evaluated to be 
negligible, these energy bins participate in the day/night 
variation. To combine the time variation constraints with 
those from the recoil electron spectrum, A log C is inter- 
preted as a time- variation Ax^^ = — 2A log C and added 
to the spectrum x^- 

Disregarding oscillations, we calculate the interaction 
rates in bin i ^B^ (hepj due to *B (hep) neutrinos. We 
also compute the oscillated rates ^B"'" and hepf for 
each tan^ 6 and Am^ . From these and the measured 
rates Data^, we form the ratios di =T)atai / {^Bi+hepi) , 
h =«Br7(8B,-fhep,), and/i, =8hepr/(''B,+hep,). The 
expected oscillation suppressions are (3bi + rjhi where 
f3 [t]) is the ^B (hep) neutrino flux scaling parameter. 
These suppressions are modified by the correlated un- 
certainty distortion functions ffiSs) (uncertainty in the 
®B neutrino spectrum), ff (Ss) (uncertainty in SK energy 



scale), and /,/^((5i?) (uncertainty in SK energy resolution) 
to p, = §h+Ilhi. where /, = fff.'^ff- With the total 
energy bin-uncorrelated uncertainty af = crf^t^t + ffay^u 
the total x^ is then 
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where the last term constraining the '^B flux to the stan- 
dard solar model (SSM) Q is optional. Including this last 
term, the best oscillation fit is in the quasi- vacuum region 
at Arn^ = 6.49 x 10^^°eV^ and maximal mixing, where 
a summer/winter asymmetry of -0.6% is expected and 
-0.3±0.7%(stat) is found. The x^ is 17.1 for 20 degrees 
of freedom (65% C.L.). The LMA solution fits almost 
equally well: the smallest x^ at Am^ = 6.3 x 10~^eV^ 
and tw^e = 0.55 is 17.3 (63% C.L.). Figure El shows 
the allowed areas at 95% C.L. using all SK information: 
rate, spectrum and time- variation. It also shows the Ax^ 
as a function of tan^ 6* (Arn^) alone: SK data excludes 
small mixing at more than 3cr. SK data also disfavors 
Am2 > lO-^eV^ and2xlO-^eV2 < Am^ < 2,xlQ-^eY'^. 

Stronger constraints on Am^ result from the combi- 
nation of SK with other solar neutrino data The 
combined fit to SK data and the SNO measurements on 
the charged-current and neutral-current reactions of so- 
lar ®B neutrinos with deuterons need not constrain any 
neutrino flux with a solar model. Figure shows the 
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FIG. 7: Ax^ Curves as a Function of Mixing (Top) and Am^ 
(Bottom) Using all Solar and KamLAND Data. Only LMA-I 
remains allowed at 3a. 



allowed region at 95%C.L.: only LMA solutions survive. 
When the charged-current rates measured by Homestake, 
GALLEX/GNO, and SAGE are included as weU, the 
LMA solutions are favored by 3a; however, the fit re- 
lies on the SSM predictions of the pp, pep, CNO, and 
^Be neutrino fluxes. 

The first oscillation analysis of the KamLAND reactor 



neutrino spectrum and rate leaves several allowed areas, 
usually called LMA-0, LMA-I, LMA-II, and LMA-III. 
When we combine the analysis of SK with all other solar 
experiments and a likelihood analysis of the KamLAND 
data 0, the LMA-I is strongly favored over the other 
solutions: in Figured the Ax^ of the fit is plotted against 
Am^ and tan^ 9. 

In summary, SK has measured very precisely the 
neutrino flux time-variations expected from two- 
neutrino oscillations. For the best LMA parameters, 
the day/night asymmetry is determined as = —1.8 ± 
1.6(stat)+};2(syst)% where -2.1% is expected. SK data 
disfavor large Arn^ LMA solutions, since their expected 
day /night asymmetries are closer to zero. In combination 
with other solar data and the KamLAND reactor neu- 
trino results, the oscillation parameters are determined 
as Am^ = 7.1+_oi x IQ-^eV^ and tan^ 9 = 0.44 ± 0.08. 
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